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ABSTRACT 



C^ We present colour-colour diagrams of detected sources in the Herschel-AiTLAS Science Demonstration Field from 100 to 500 yum using both 

r ■) PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with 

^*^ counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 

r^ 330 sources detected in at least three FIR bands with a significance greater than 3o", we find an average dust temperature of (28 ± 8)K. For sources 

O I with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures 

and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample 

of 1686 sources with fluxes above 35 mjy at 350;/m and detected at 250 and 500 ;/m with a significance greater than 3cr, we find an average 

redshift of 2.2 ± 0.6. 

C^ Key words, submillimeter: galaxies — Galaxies: evolution — Galaxies: high-redshift 
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>. ■ 1 . Introduction ciding with the GAMA redshift survey fields (Driver et al. 2009), 

CNj ' and 2 fields (total area of 250 deg^) near the southern Galactic 

,__l The total intensity of the extragalactic background light at far- pQjg 

^ infrared (FIR) wavelengths has been measured with absolute ^^^.- ^^^ ^^-^^^^ Demonstration Phase (SDP), H- ATLAS 

(N photometry (Puget et al. 1996; Fixsen et al. 1998; Dwek et al. observed 14.4 deg^ in the GAMA-9 hour field near the ecliptic 

^^ 1998). Deep surveys are now starting to resolve the cosmic FIR pj^^^^ ^^ ^ 5^ ^^p^j^^ ^f 35.9O ^j^ j^ ^^is Letter we present the 

j-5 background into discrete sources with the fraction resolved vary- colour-colour diagrams derived from the five bands from 100 to 

O ing with wavelength (see reviews in Blain et al. 2002; Hauser & ^qq^^ ^^^ jj^^ H-ATLAS sources that are cross-identified in at 

^ , Dwek 2001 ; Lagache et al. 2005, Dye et al. 2009). Despite some j^^^^ 3 ^f jj^^ ^pj^g ^^^ p^^S bands. We then discuss the dust 

LJ ■ successes, there are still large uncertainties about the nature and temperature and redshift distributions based on certain assump- 

evolution of the submillimetre (submm) galaxy population, in- ^-^^^ ^^^^^ ^j^^ ^p^^^^^j ^^^^gy distribution (SED) and subject to 

eluding their redshift distribution. selection eff^ects associated with the sample selection. 
^ ] The i/er5c/ieZ-AstrophysicaI Terahertz Large Area Survey 

(H-ATLAS) is an open-time key program of the Herschel 

Observatory (Pilbratt et al. 2010) that will survey roughly 550 

deg^ of sky over 600 hours of observations (Eales et al. 2010). 
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2. Colour-colour diagrams 



To cover the widest possible area, H-ATLAS uses the maximum We use source catalogues generated for the H-ATLAS consor- 

possible scan rate for the telescope at 60 arcsec sec ' in the par- tium by Rigby et al. (2010, in prep.), derived from SPIRE and 

allel mode of PACS (Poglitsch et al. 2010) and SPIRE (Griffin PACS maps presented in Pascale et al. (2010, in prep.) and 

et al. 2010). The observations cover five photometric bands with Ibar et al. (2010, in prep.), respectively. The source catalogues 

100 and 160 yum data from PACS, and 250, 350, and 500 yum data are supplemented with cross-identification information from the 

from SPIRE. The H-ATLAS fields include 1 field close to the GAMA survey (Driver et al. 2009) and SDSS DR-7 (Abazajian 

northern Galactic pole (150 deg^), 3 fields (each 36 deg^) coin- et al. 2009) as described in Smith et al. (2010, in prep.). For 

the sources for which spectroscopic redshifts from GAMA or 

* Herschel is an ESA space observatory with science instruments SDSS are not available, we use photometric redshifts gener- 

provided by European-led Principal Investigator consortia and with im- ated using the ANNz neural network code (Collister & Lahav 

portant participation from NASA. 2004), trained with photometry from SDSS and UKIDSS LAS 
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Fig. 1. Colour-colour diagram of SPIRE and PACS sources. 
The coloured background indicates the average redshift in these 
colour-colour spaces of our 10^ randomly generated model 
SEDs. In (a) points plotted represent a sample of 1686 galax- 
ies detected in all three SPIRE bands with a significance greater 
than 5 cr at 350 ;um (S > 35mJy) and greater than 3cr at 250 and 
500 //m (coiTesponding to flux limits of roughly 21 and 27 mJy, 
respectively). In (b) we show a sample of 402 galaxies detected 
at 160, 250, and 350/vm and in (c) we show a sample of 158 
galaxies detected at 100, 160, and 250 yum. As in (a), in both (b) 
and (c) we also select sources by imposing a 5cr cut in the cen- 
tral wavelength of the colour plot and greater than 3cr at other 
wavelengths. 



(Lawrence et al. 2007), and spectra from the GAMA spectro- 
scopic survey (Driver et al 2009), DEEP2 (Davis et al. 2007), 
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Fig. 2. Dust temperature as a function of redshift for 330 H- 
ATLAS sources reliably identified with GAMA and SDSS DR- 
7 galaxies having a known spectroscopic (black filled triangle) 
or a photometric (black empty triangle) redshift. The selected 
sources are > 5cr detections in at least one band of either PACS 
or SPIRE with 3cr flux measurements in at least 2 more bands. 
The isothermal spectral fits assume/? = 1.5. For other data plot- 
ted, see Section 3. 
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Fig. 3. Dust temperature as a function of the FIR luminosity 
for 330 H-ATLAS sources reliably identified with GAMA and 
SDSS DR-7 galaxies with a known spectroscopic (black filled 
triangle) or a photometric (black empty triangle) redshift. The 
FIR luminosity is obtained by fitting the PACS and SPIRE mea- 
surements with a modified black body (parameters are the dust 
temperature and the luminosity density,/? is fixed to 1.5) and this 
model SED is integrated between 8 and 1 lOOyum. We also indi- 
cate BLAST and H-ATLAS best-fit Td - \og(LFiR) relation with 
a dashed blue line (from Dye et al. 2009) and a dashed red line, 
respectively. 



zCOSMOS (Lilly et al. 2007), and the 2SLAQ-LRG (Cannon et 
al. 2006) survey. 

In Fig. 1, we show the colour-colour plots of the H-ATLAS 
sources. We divide these plots in terms of colours based on either 
SPIRE only or SPIRE and PACS data. Our flux selection results 
in selecting 1686, 402, and 158 sources, from top to bottom of 
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Fig. 1. For reference, the total ATLAS catalogue for this field 
contains ~ 6600 sources (Rigby et al. 2010, in prep.). 

The colour-colour plots are filled with lO*" black-body spec- 
tra at a single dust temperature, T^, modified by a frequency- 
dependent emissivity function €y oc v", where the flux density fy 
is 



2.0 



/v = evBv'^v3+/'/[exp(^)-l]. 



(1) 



In generating these models, we consider uniform ranges of dust 
temperature from 10 to 60 K, emissivity parameter's from to 

2, and redshift from to 5. The choice of < yS < 2 compared to 
1 < j0 < 2 makes a minor difference since we also broaden the 
SED tracks in the colour diagram by adding an extra Gaussian 
standard deviation of 10% to the fluxes used to compute the 
model colour This scatter accounts for the broadening of data in 
the colour-colour plane caused by flux uncertainties. As shown 
in Fig. 1, we find that the colour diagram of sources with fluxes 
only from SPIRE are well within the limits defined by the mod- 
els we have considered. When we examine PACS colours, we 
find that some points lie outside the same set of tracks as used 
for the SPIRE-only colour diagram. While some of these outlier 
points may be cause by either the fractionally larger flux errors 
of PACS or contamination from a neighboring source, it is pos- 
sible that some of these sources are not accurately described by 
our simple isothermal SED model, requiring for instance a sec- 
ond dust component (Dunne & Eales 2001) or a more complex 
SED model. 

When fitting a simple modified black-body model to the data, 
we must keep in mind that there is a partial degeneracy between 
y6 and Tj and, more importantly, a perfect degeneracy between 
Td and z- The peak of the SED is determined by the v/T^ term in 
the exponential, so that a measurement of the colours alone con- 
strains only the ratio (I + zj/T^. However, assuming reasonable 
priors on the free parameters of the SED model (in our case, 
j3 and Tj) it is still possible to estimate a qualitative redshift 
distribution for our sample of sources. Alternatively, if secure 
redshifts are known from optical cross-identifications, we can 
determine the dust temperatures. 

3. Dust temperature distribution 

We consider the H-ATLAS source sample with detections at 3cr 
in at least two bands and at 5cr in one band of either PACS and 
SPIRE that have been robustly (reliability parameter Rlr > 0.9, 
Smith et al. 2010, in prep.) identified with GAMA or SDSS DR- 
7 galaxies. We also require that there is a known spectroscopic 
redshift from either GAMA or SDSS, or from the photometric 
redshift catalog (with (1 + z)/o-- > 5 and z/cr- > 1) that was 
generated for this field and cross-identified with ATLAS sources 
(Smith et al. 2010, in prep.). We select 330 sources, which cor- 
respond to the low redshift subsample of the galaxies selected in 
the colour-colour diagrams. For each galaxy, we perform a sin- 
gle temperature fit from the above equation assuming yS = 1.5. 
We found that an isothermal SED model generally is a good fit to 
the 330 galaxies. Fitting a two-component SED model does not, 
on average, seem to provide a closer fit. The relatively good fit 
provided by the isothermal SED model is probably due to most 
of our sample consisting of low redshift galaxies and therefore 
that we do not probe much of the Wien part of the galaxy SED. 
In Fig. 2, we summarize our results and compare the H- 
ATLAS dust temperatures with dust temperatures in the liter- 
ature for a variety of sub-mm bright galaxies. These samples 
are: (i) the sources in BLAST detected above 5cr in at least 
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Fig. 4. Normalized redshift distribution (black thick solid line) 
of SPIRE sources in the colour-colour diagram of Fig. 1(a), se- 
lected to satisfy the following criteria: S 350 > 35mJy and at least 
3cr detections at 250 and 500 yum. The dashed line and dotted- 
dash line represent Negrello et al. (2007) and Lagache et al. 
(2004) models for the same flux selection. Both model curves 
have similar but not identical distributions, although Lagache 
model predicts a larger number of high-redshift galaxies. In ad- 
dition, we show N{z) for the same subsample using a colour-cut 
of 5500/5 250 < 0.75 (blue thin solid line), and > 0.75 (red thin 
dotted line). The distribution shown here is for 1686 sources, a 
subset of about 6600 H-ATLAS sources in the SDP field. 



one of the BLAST bands with either a COMBO-17 (Wolf et al. 
2004) or a SWIRE photometric redshift (Rowan-Robinson et al. 
2008) and Spitzer-MlPS 70 and 160 //m fluxes (Dye et al. 2009; 
yS = 1.5 fixed); (ii) local ULIRGS observed with SCUBA at 450 
and 850 //m and complemented with IRAS 60 and lOOyum fluxes 
(Clements et al. 2010; y6 varied); (iii) SCUBA sub-mm galaxies 
detected at better than 3cr at 850 /im and having redshifts deter- 
mined from Keck-I spectroscopy (Chapman et al. 2005; y6 =1.5 
fixed); and (iv) local IRAS-selected galaxies with 60 and 100 /im 
fluxes complemented with SCUBA 850 yum (Dunne et al. 2000; 
j8 varied). 

In Table 1, we list average dust temperatures as a function 
of redshift for several bins in redshift for both H-ATLAS only 
and all of the combined sub-mm galaxy samples, including H- 
ATLAS, plotted in Fig. 2. We do not find any evolution in the H- 
ATLAS dust temperature with redshift, though some evolution 
is inferred by BLAST measurements (Dye et al. 2009; Pascale 
et al. 2009), where sources at higher redshift had higher tem- 
perature in agreement with the radio-identified submillimeter- 
selected galaxies (SMGs) found with SCUBA (Chapman et al. 
2005; Ivison et al. 2010; Kovacs et al. 2006; Coppin et al. 2008). 
While we find remarkable consistency between the average dust 
temperatures for our H-ATLAS sample as a function of redshift, 
these average values are not necessarily in agreement with other 
sub-mm galaxy subsamples in the literature, mostly due to se- 
lection eff'ects. For example, the SCUBA ULIRGS have an av- 
erage temperature of (43 + 7)K. These sources were selected us- 
ing IRAS 6Qfim data and IRAS selected sources are known to 
be biased towards higher temperatures. Optically-selected low- 
redshift sub-mm galaxies are known to have colder temperatures 
consistent with our findings (e.g., Willmer et al. 2009; Vlahakis 
et al. 2005). The large expected sample of sources from the 550 
deg.^ of H-ATLAS will enable more detailed studies in the fu- 
ture, without the biases associated with selections of various sub- 
samples, including our own. While we show results here with 
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Table 1. Average dust temperatures as a function of redshift for 
the 330 H-ATLAS galaxies (Col. 2) and for all the data (Col. 4) 
presented in Fig. 2 and 3 (including H-ATLAS). 



z-range 


H-ATLAS 

Ns,-cs 


Td 


all data 

Ns,cs 


all data T,] 


Allz 


330 


28 ±8 


657 


30 ±9 


0<z<0.1 


106 


27 ±8 


235 


32 ±9 


0.1 <z<0.5 


186 


29 ±8 


260 


28 ±8 


0.5 < z < 1 


33 


23 ±5 


67 


24 ±7 


z> 1 


5 


35 ±4 


95 


37 ±10 



j3 - 1.5 fixed, when fitting foryS and T^ we found y8 = 1.4 + 0.1, 
consistent withyS =1.3 found in Dunne et al. (2000). 

In Fig. 3, we plot the dust temperature versus FIR lumi- 
nosity by integrating model SEDs between 8 and 1100 fim for 
the H-ATLAS subsample. Luminosities for other samples are 
from the literature. Fitting for a relation of the form T^ - 
To + alogiLfiRlLo), we find Tq - -20.5K and a - 4.4 (see. 
Fig. 4). The value oflogiLfiRlLo) is on average 10.9+0.8 for our 
sample. This relation is consistent with the BLAST data (Dye et 
al. 2009). 



4. Redshift distribution 

Since most of our sources lack redshifts, we consider another 
example application of our colour diagram and infer the statis- 
tical redshift distribution N(z) for the source samples plotted in 
Fig. 1. We do this by first gridding the colour-colour plane along 
the SED tracks into redshift bins. We then convert the number 
of sources within a grid region of colours to a binned redshift 
distribution (Hughes et al. 2002). This method is equivalent to 
extracting the redshift probability distribution function for the 
whole sample if we had simply fitted SEDs to individual fluxes 
and taken the sum of the redshift probabilities of each source. 
While the redshift for an individual galaxy is largely uncertain, 
and sensitive to the SEDs used, the statistical redshift distribu- 
tion we extract should be a reasonable estimate of the true dis- 
tribution of the source sample. 

Figure 4 shows the redshift distribution for SPIRE sources 
detected in all 3 bands (Fig. 1(a)). For the sample of 1686 
sources with flux densities above 35 mJy at 350 /zm and above 
3cr at 250 and 500 ;um, we find the average redshift to be 2.2+0.6. 
This is consistent with the average redshift of 2.4 + 0.4 for the 
radio-identified sub-mm galaxies with SCUBA at 850 jum fol- 
lowed up spectroscopically at Keck (Chapman et al. 2005), the 
interquartile redshift range of 1.8 to 3.1 for sources in SCUBA 
Half Degree Extragalactic Survey determined photometrically 
(Aretxaga et al. 2007), and a median redshift of greater than 1 
for a subsample of 250 yum sources as faint as 35 mJy (Dunlop 
et al. 2009). Our distribution, however, has an average redshift 
that is higher than that of a larger sample of BLAST sources 
(Dye et al. 2009). This difference is not unexpected given the 
differences in the sample selection. We emphasize again that 
the distribution shown in Fig. 4 is for a subset of about 25% 
of the whole catalog with specific selection criteria that tend to 
favour the high redshift-end of the distribution. Based on cross- 
identifications with SDSS and GAMA, we find that 30% of the 
H-ATLAS sources are at z < 0.5 in the SDP field (Smith et al. 
2010, in prep.). 

In addition to the total sample, we also make a colour-cut 
with S500/S250 > 0.75 and < 0.75, in addition to the flux cut of 
35 mJy at 350;um. By imposing this colour-cut, we distinguish a 



high redshift population with redder colours from a low redshift 
population. We find average redshifts of 2.6+0.3 and 1.8+0.4 for 
sources with S;,oo/S250 > 0.75 and < 0.75, respectively. While 
we do not show the redshift distributions separately for the sam- 
ples in Figs. 1(b) and (c), the average redshift in these cases is 
z = 0.8 + 0.3 for sources detected in 160, 250, and 350;[/m with 
a flux cut of 35 mJy at 250 /im, and z = 0.7 + 0.3 for sources 
detected in 100, 160, and 250 //m with a flux cut of ^ 100 mJy 
at 160 /zm. As is clear from Fig. 1, we find that lower redshift 
sources are more likely to be detected at shorter wavelengths, 
while bright sources seen in longer SPIRE bands are likely to be 
on average at high redshifts. Some of the bright 500 /im sources 
that we detect, especially with a colour selection selecting red 
galaxies in the sub-mm, may also have their fluxes magnified 
by lensing and will be among the results forthcoming from H- 
ATLAS (Negrello et al. 2010, in prep.). 

5. Conclusions 

We have discussed the spectral energy distribution of sub-mm 
galaxies in a 14 deg.^, subset of the full H-ATLAS survey. We 
have used the colours in 5 bands to measure the dust temper- 
ature for the sub-sample of sources with known redshifts. We 
have also derived a qualitative estimate of the statistical redshift 
distribution of bright 350 /im selected galaxies assuming broad 
priors in temperature and the dust emissivity. We found that the 
average dust temperature of H-ATLAS sources is 28 ± 8K, while 
350 /im sources with fluxes above 35 mJy have an average red- 
shift of 2.2 + 0.6. 
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